Effective temperature and Gilbert damping of a current-driven localized spin 
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Starting from a model that consists of a semiclassical spin coupled to two leads we present a 
microscopic derivation of the Langevin equation for the direction of the spin. For slowly-changing 
direction it takes on the form of the stochastic Landau-Lifschitz-Gilbert equation. We give ex- 
pressions for the Gilbert damping parameter and the strength of the fluctuations, including their 
bias-voltage dependence. At nonzero bias-voltage the fluctuations and damping are not related by 
the fluctuation-dissipation theorem. We find, however, that in the low-frequency limit it is possible 
to introduce a voltage-dependent effective temperature that characterizes the fluctuations in the 
direction of the spin, and its transport-steady-state probability distribution function. 
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I. INTRODUCTION 

One of the major challenges in the theoretical descrip- 
tion of various spintronics phenomena 1 , such as current- 
induced magnetization reversal 2 ^^ 5 - and domain-wall 
motion 6 -! 7 ^ ' 9 ' 10 ' 11 ' 12 , is their inherent nonequilibrium 
character. In addition to the dynamics of the collective 
degree of freedom, the magnetization, the nonequilibrium 
behavior manifests itself in the quasi-particle degrees of 
freedom that are driven out of equilibrium by the nonzero 
bias voltage. Due to this, the fluctuation-dissipation 
theore m 13 ' 14 cannot be applied to the quasi-particles. 
This, in part, has led to controversy surrounding the the- 
ory of current-induced domain wall motio n 15 : 16 . 

Effective equations of motion for order-parameter 
dynamics that do obey the equilibrium fluctuation- 
dissipation theorem often take the form of Langevin 
equations, or their corresponding Fokker-Planck 
equation a 13 ' 14 : 17 . In the context of spintronics the rele- 
vant equation is the stochastic Landau-Lifschitz-Gilbert 
equation for the magnetization directio n 18 : 19 ' 20 : 21 : 22 : 23 : 24 . 
In this paper we derive the generalization of this equa- 
tion to the nonzero-current situation, for a simple 
microscopic model consisting of a single spin coupled 
to two leads via an onsite Kondo coupling. This model 
is intended as a toy-model for a magnetic impurity 
in a tunnel junctio n 25 ' 26 ' 27 . Alternatively, one may 
think of a nanomagnet consisting of a collection of 
spins that are locked by strong exchange coupling. 
The use of this simple model is primarily motivated 
by the fact that it enables us to obtain analytical 
results. Because the microscopic starting point for 
discussing more realistic situations has a similar form, 
however, we believe that our main results apply quali- 
tatively to more complicated situations as well. Similar 
models have been used previously to explicitly study 
the violation of the fluctuation-dissipation relation 28 , 
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FIG. 1: Effective temperature as a function of bias voltage. 
The dashed line shows the large bias-voltage asymptotic result 
teTeff — |e| V/4 + kBT/2. The inset shows the bias-voltage 
dependence of the Gilbert damping parameter normalized to 
the zero-bias result. 



and the voltage-dependence of the Gilbert damping 
parameter—. Starting from this model, we derive an 
effective stochastic equation for the dynamics of the 
spin direction using the functional-integral description 
of the Keldysh-Kadanoff-Baym nonequilibrium theory 29 . 
(For similar approaches to spin and magnetization 
dynamics, see also the work by Rebei and Simionatc 30 , 
Nussinov et aL— and Duine et o/.— .) This formalism 
leads in a natural way to the path-integral formulation 
of stochastic differential equation a 33 ' 34 . One of the 
attractive features of this formalism is that dissipation 
and fluctuations enter the theory separately. This allows 
us to calculate the strength of the fluctuations even 
when the fluctuation-dissipation theorem is not valid. 

We find that the dynamics of the direction of the spin 
is described by a Langevin equation with a damping ker- 
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FIG. 2: Model system of a spin S connected to two tight- 
binding model half-infinite leads. The chemical potential of 
the left lead is /xl and different from the chemical potential 
of the right lead /^r. The temperature T of both leads is for 
simplicity taken to be equal. 



nel and a stochastic magnetic field. We give explicit 
expressions for the damping kernel and the correlation 
function of the stochastic magnetic field that are valid 
in the entire frequency domain. In general, they are not 
related by the fluctuation-dissipation theorem. In the 
low-frequency limit the Langevin equation takes on the 
form of the stochastic Landau-Lifschitz-Gilbert equation. 
Moreover, in that limit it is always possible to introduce 
an effective temperature that characterizes the fluctua- 
tions and the equilibrium probability distribution for the 
spin direction. In Fig. [1] we present our main results, 
namely the bias- voltage dependence of the effective tem- 
perature and the Gilbert damping parameter. We find 
that the Gilbert damping constant initially varies lin- 
early with the bias voltage, in agreement with the re- 
sult of Katsura et al£L. The voltage-dependence of the 
Gilbert damping parameter is determined by the den- 
sity of states evaluated at an energy equal to the sum 
of the Fermi energy and the bias voltage. The effective 
temperature is for small bias voltage equal to the actual 
temperature, whereas for large bias voltage it is inde- 
pendent of the temperature and proportional to the bias 
voltage. This bias-dependence of the effective tempera- 
ture is traced back to shot noised. 

Effective temperatures for magnetization dynam- 
ics have been introduced before on phcnomcnolog- 
ical grounds in the context of thermally-assisted 
current-driven magnetization reversal in magnetic 
nanopillar o 36 ' 37 i 38 . A current-dependent effective tem- 
perature enters in the theoretical description of these 
systems because the current effectively lowers the energy 
barrier thermal fluctuations have to overcome. In addi- 
tion to this effect, the presence of nonzero current alters 
the magnetization noise due to spin current shot noise 3 -^. 
Covington et al. 39 interpret their experiment in terms of 
current-dependent noise although this interpretation is 
still under debate 3 - . Foros et al£5. also predict, using a 
different model and different methods, a crossover from 
thermal to shot-noise dominated magnetization noise for 
increasing bias voltage. Our main result in Fig. [1] is an 
explicit example of this crossover for a specific model. 

The remainder of the paper is organized as follows. We 
start in Sec.|TT]by deriving the general Langevin equation 
for the dynamics of the magnetic impurity coupled to two 
leads. In Sec. IIIII and IIVI we discuss the low-frequency 
limit in the absence and presence of a current, respec- 
tively. We end in Sec. [V] with our conclusions. 



II. DERIVATION OF THE LANGEVIN 
EQUATION 

We use a model that consists of a spin S on a site 
that is coupled via hopping to two semi-infinite leads, as 
shown in Fig. [2] The full probability distribution for the 
direction Cl of the spin on the unit sphere is written as a 
coherent-state path integral over all electron Grassmann 
field evolutions tp*(t) and tp(t), and unit-sphere paths 
S(t), that evolve from — oo to t and back on the so-called 
Keldysh contour C* . It is given by22i 



P[Q,t] 



d[S] S 

s(t)=n 
xexp(-^*,^,S] 



(1) 



where the delta functional enforces the length constraint 
of the spin. In the above functional integral an inte- 
gration over boundary conditions at t = — oo, weighted 
by an appropriate initial density matrix, is implicitly in- 
cluded in the measure. We have not included boundary 
conditions on the electron fields, because, as we shall see, 
the electron correlation functions that enter the theory 
after integrating out the electrons are in practice conve- 
niently determined assuming that the electrons are cither 
in equilibrium or in the transport steady state. 
The action S[ip*,ip, S] is the sum of four parts, 
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(2) 



We describe the leads using one-dimensional non- 
interacting electron tight-binding models with the action 



L/R 



^lE(e(*'))*4e(o 



3,a 



(j,j'h 



(3) 



where the sum in the second term of this action is 
over nearest neighbors only and proportional to the 
nearest- neighbor hopping amplitude J in the two leads. 
(Throughout this paper the electron spin indices are de- 
noted by (J, (j 1 e {t; 4}, and the site indices by j, j'.) The 
coupling between system and leads is determined by the 
action 



3 



(4) 



where dL and <9R denote the end sites of the semi-infinite 
left and right lead, and the fields (ip°(t)) and ip°(t) de- 
scribe the electrons in the single-site system. The hop- 
ping amplitude between the single-site system and the 
leads is denoted by Jo Finally, the action for the sys- 
tem reads 



S° 



= / dt' 

Jc* 



d 



dt' 



A£(v£(*oy 



s(f) rA* 



(5) 



The second term in this action is the usual Berry phase 
for spin quantization^, with A(S) the vector potential 
of a magnetic monopole 



ds, 
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(6) 



where a sum over repeated Greek indices a,(3,j <E 
{x,y,z} is implied throughout the paper, and t a p-y is 
the anti-symmetric Levi-Civita tensor. The third term in 
the action in Eq. ([5|) describes the coupling of the spin to 
an external magnetic field, up to dimensionful prefactors 
given by h. (Note that h has the dimensions of energy.) 
The last term in the action models the s — d exchange 
coupling of the spin with the spin of the conduction elec- 
trons in the single-site system and is proportional to the 
exchange coupling constant A > 0. The spin of the con- 
duction electrons is represented by the vector of the Pauli 
matrices that is denoted by t. 

Next, we proceed to integrate out the electrons using 
second-order perturbation theory in A. This results in 
an effective action for the spin given by 



S ctt [S] 



dt' 



ShA(S(t')) 



dS(t') 



dt' 



■h-S(i') 



dt"U(t',t")S(t') ■ S(t") 



(7) 



This perturbation theory is valid as long as the electron 
band width is much larger than the exchange interac- 
tion with the spin, i.e., J, Jq 3> A. The Keldysh quasi- 
particle response function is given in terms of the Keldysh 
Green's functions by 



-~G(t,t')G(t',t) 



n(M') 

where the Keldysh Green's function is defined by 



(8) 



(9) 



We will give explicit expressions for the response function 
and the Green's function later on. For now, we will only 
make use of the fact that a general function A(t, t') with 
its arguments on the Keldysh contour is decomposed into 
its analytic pieces by means of 



A(t, t') = 9(t, t')A> (t, t') + 6(t', t)A< (t, t') 



(10) 



where 9(t, t') is the Heaviside step function on the 
Keldysh contour. There can be also a singular piece 
A s S(t, t') 1 but such a general decomposition is not needed 
here. Also needed are the advanced and retarded com- 
ponents, denoted respectively by the superscript (— ) and 
(+), and defined by 

A^{t,t') = ±e{±{t-t'))[A>{t,t')-A<{tJ)\ , (11) 

and, finally, the Keldysh component 

A K (t,t')=A > (t,t')+A < (t,t'), (12) 

which, as we shall see, determines the strength of the 
fluctuations. 

Next we write the forward and backward paths of the 
spin on the Keldysh contour, denoted respectively by 
S(£ + ) and S(£_), as a classical path Cl(t) plus fluctua- 
tions Stt(t), by means of 



s(f±) = n(t)± 



sn(t) 



(13) 



Moreover, it turns out to be convenient to write the delta 
functional, which implements the length constraint of the 
spin, as a path integral over a Lagrange multiplier A(t) 
defined on the Keldysh contour. Hence we have for the 
probability distribution in first instance that 



p[Ci,t] = 



s(t)=o 



with 



d[SMA]exp<j-S ctt [S] + -S A [S,A]j , 

(14) 



5 A [S,A] = / dt'A(t') [|S(0| 2 -1 



(15) 



We then also have to split the Lagrange multiplier into 
classical and fluctuating parts according to 



A(t±) = A(t) ± 



5X(t) 



(16) 



Note that the coordinate transformations in 
Eqs. (fl"3"|) and (fTI)]) have a Jacobian of one. Before 
we proceed, we note that in principle we are required 
to expand the action up to all orders in Stt. Also note 
that for some forward and backward paths S(£+) and 
S(£_) on the unit sphere the classical path S~2 is not 
necessarily on the unit sphere. In order to circumvent 
these problems we note that the Berry phase term in 
Eq. ([5]) is proportional to the area on the unit sphere 
enclosed by the forward and backward paths. Hence, in 
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the semi-classical limit S — ► oo 27 i 40 paths whose forward 
and backward components differ substantially will be 
suppressed in the path integral. Therefore, we take this 
limit from now on which allows us to expand the action 
in terms of fluctuations 6fl(t) up to quadratic order. We 
will see that the classical path fl(t) is now on the unit 
sphere. We note that this semi-classical approximation 
is not related to the second-order perturbation theory 
used to derive the effective action. 

Splitting the paths in classical and fluctuation parts 
gives for the probability distribution 



- A 2 
f 



dt / dt 



t'{sn a (t) [n(-\t',i)+n(+\M')] fi a (*0} 



2A 2 



dt / dt' 



v a (t)(n K ) 1 (t,t') Va (t>) 



(21) 



Note that the inverse in the last term is defined as 

/ dt"u K (t : t") (n K ) _1 (t", f ) = 8(t - if). 

Performing now the path integral over 6fl(t), we ob- 
serve that the spin direction fl(t) is constraint to obey 
the Langevin equation 



hSe a (3 7 — ^— -fi 7 (i) 



p[ti,t] = 



d[n}d[5n}d[x}d[sx\ ex P <^ -s[n, sti, x, sx\ 



(17) 

with the action, that is now projected on the real-time 
axis, 



h a + 2X(t)fl a (t) 

OO 

dt'K(t,t')n a (t') , (22) 



with the so-called damping or friction kernel given by 



S[ft, Sfl, A, SX] = J dt IhSeafoSnpfy) 
+5Vt a (t)h a + 2SCl a (t)Cl a (t)X(t) 



x dn a (t) 

~dt 



K(t,t') 



n(-)(t',t) + n( +) (t,i') 



(23) 



n 7 (t) 



-5X{t) [\fl(t)\ 



\sn(t)\ 2 /A] 



Note that the Heaviside step functions in Eq. (fTTjl appear 
precisely such that the Langevin equation is causal. The 
stochastic magnetic field is seen from Eq. (l2TT) to have 
the correlations 



dt / dt' 



n(- 



A 2 

A 2 

— I dt I dt' [sn a (t)u K (t,t')sn a (t')] 



(v a (t)) = o ; 

>(t,0] n a {t')} (va(t)v0(t')) = t<y a /sfiA 2 n K (t,0 . (24) 

Using the fact that H(t) is a unit vector within our semi- 
(18) classical approximation, the Langevin equation for the 
direction of the spin Sl(t) is written as 



From this action we observe that the integration over 
SX(t) immediately leads to the constraint 



hS 



dfi(t) 

dt 



n{t) x 



h + 77(i) 



dt'K(t,t')Cl{t') 



\n(t)\ 



1 



\6Sl(t)\< 



(19) 



as expected. Implementing this constraint leads to terms 
of order 0(5tt 3 ) or higher in the above action which we 
are allowed to neglect because of the semi-classical limit. 
From now on we can therefore take the path integration 
over fi(t) on the unit sphere. 

The physical meaning of the terms linear and quadratic 
in 6Sl(t) becomes clear after a so-called Hubbard- 
Stratonovich transformation which amounts to rewrit- 
ing the action that is quadratic in the fluctuations as 
a path integral over an auxiliary field rj{t). Performing 
this transformation leads to 



(25) 

which has the form of a Landau-Lifschitz equation with 
a stochastic magnetic field and a damping kernel. In 
the next sections we will see that for slowly- varying spin 
direction we get the usual form of the Gilbert damping 
term. 

So far, we have not given explicit expressions for the 
response functions n( ± ' ,K (t, t'). To determine these func- 
tions, we assume that the left and right leads are in 
thermal equilibrium at chemical potentials /il and /zr, 
respectively. Although not necessary for our theoretical 
approach we assume, for simplicity, that the temperature 
T of the two leads is the same. The Green's functions for 
the system are then given b y 41 ' 42 



p[Ci,t] = 



d[Ci\d[5Sl]d[rj\d[X] 
xexp\^S[fl,5fl,X 7 r]} 



(20) 



where the path integration over Q is now on the unit 
sphere. The action that weighs these paths is given by 



-iG<(e) = 

ZG>(6) = 

G^' K (t-t') = 



2 

fce{L,R} 



s[n,«i,A,»j] = J 



2 E 

fce{L,R} 
de 



[l-N(e-fi k )] 



dt 



hSe af}7 SQ p (t)^^Q y (t) 



(27T) 



e -u(f-tfyhQ$^ ^ (26) 



with N(e) = {exp[e/(fc B 7 1 )] + the Fermi-Dirac dis- 
tribution function with fee Boltzmann's constant, and 



-5n a (t)h a + 2sn a (t)n a (t)x(t) + sn a (t)r] a (t) 



A(e) =1 [G (+) (e)-G ( - ) (e) 



(27) 
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the spectral function. Note that Eq. ([26)1 has a particu- 
larly simple form because we are dealing with a single-site 
system. The retarded and advanced Green's functions 
are determined by 



e ± -2h^ ± \e)] G (±) (e) = 1 



(28) 



with e^ 1 = e ± iO, and the retarded self-energy due to one 
lead follows, for a one-dimensional tight-binding model, 



J, 



/i£ (+) (e) = —j-e 



C_ ik(e)a 



(29) 



with k(e) = arccos[— e/(2 J)]/a the wave vector in the 
leads at energy e, and a the lattice constant. The ad- 
vanced self-energy due to one lead is given by the com- 
plex conjugate of the retarded one. 

Before proceeding we give a brief physical description 
of the above results. (More details can be found in 
Refs. [4l[ and [Hj].) They arise by adiabatically elim- 
inating ("integrating out") the leads from the system, 
assuming that they are in equilibrium at their respective 
chemical potentials. This procedure reduces the problem 
to a single-site one, with self-energy corrections for the 
on-site electron that describe the broadening of the on- 
site spectral function from a delta function at the (bare) 
on-site energy to the spectral function in Eq. (f27f . More- 
over, the self-energy corrections also describe the non- 
equilibrium occupation of the single site via Eq. (|26p 

For the transport steady-state we have that 
n( ± )^ K (t,t') depends only on the difference of the time 
arguments. Using Eq. flSJ) and Eqs. JTDJl, (HH), and (12|) 
we find that the Fourier transforms are given by 

nW(e) = J d(t~t')e teit - t ' )/h Il (±) (t,t') 
de' [ de" 1 



and 



, (2tt) J (2tt) e± + e' - e" 
x[G < (e')G > (e")-G > {e')G < (e")] , (30) 



1 ' J (2tt) J (2tt) 



6(e + e'- e") 



x[G > (e')G < (e") + G<(e')G > (e")] . (31) 

In the next two sections we determine the spin dynamics 
in the low- frequency limit, using these expressions to- 
gether with the expressions for G>(e). We consider first 
the equilibrium case. 



III. EQUILIBRIUM SITUATION 

In equilibrium the chemical potentials of the two leads 
are equal so that we have /^l = A^R = /•*• Combining re- 
sults from the previous section, we find for the retarded 



and advanced response functions (the subscript "0" de- 
notes equilibrium quantities) that 



n^(e) 



de 1 



de 1 ' 



(2.) J (2^^ 
[N(e> - n) - N(e" fi)] 



(32) 



The Keldysh component of the response function is in 
equilibrium given by 



n£(e) 



-2m 



de' 



de" 



(2n)J (2tt) 



A(e')A(e")S(e 



{[l-iV(e'- M )] AT( e "- M )+AT(e'- M ) [l-iV(e" -/*)]} .(33) 

The imaginary part of the retarded and advanced re- 
sponse functions are related to the Keldysh component 
by means of 



n*(e) = ±2% [2JV B (e) + 1] Im U^\e) 



(34) 



with N B (e) = {exp[e/(fc B T)] - l}" 1 the Bose distribu- 
tion function. This is, in fact, the fluctuation-dissipation 
theorem which relates the dissipation, determined as we 
shall see by the imaginary part of the retarded and 
advanced components of the response function, to the 
strength of the fluctuations, determined by the Keldysh 
component. 

For low energies, corresponding to slow dynamics, we 
have that 



nW( £ )^n^(0)T^A 2 ( M ) £ 



(35) 



With this result the damping term in the Langevin equa- 
tion in Eq. ([2"5]) becomes 



dt'K(t,t')Cl(t') = - 



2tt 



dt 



(36) 



where we have not included the energy-independent part 
of Eq. (|35|) because it does not contribute to the equation 
of motion for Cl(t). In the low-energy limit the Keldysh 
component of the response function is given by 



U f(e) = ^k B T . 



(37) 



Putting all these results together we find that the dy- 
namics of the spin direction is, as long as the two leads 
are in equilibrium at the same temperature and chemical 
potential, determined by the stochastic Landau-Lifschitz- 
Gilbert equation 



HS 



dn(t) 



Q(t) X [h + r](t)] - ha n X 



d£l(t) 



dt dt 1 

with the equilibrium Gilbert damping parameter 

A 2 A 2 ((i) 



a 



2-k 



(38) 



(39) 
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Using Eqs. ([Ml), l[37|). and (39]) we find that the strength 
of the Gaussian stochastic magnetic field is determined 
by 



{■n«(t)rif)(t')) = 2a hk B TS{t - t')5 af} 



(40) 



Note that these delta-function type noise correlations 
are derived by approximating the time dependence of 
IT K (f , t') by a delta function in the difference of the time 
variables. This means that the noisy magnetic field 77(f) 
corresponds to a Stratonovich stochastic proces o 13 i 14 ' 17 . 

The stationary probability distribution function gen- 
erated by the Langevin equation in Eqs. ([381) and ([40]) is 
given by the Boltzmann distributio n 18 i 19 ' 20 i 21 ' 22 i 23 ' 24 



P[Cl, t — > 00] oc exp • 



E(Q) 
k B T 



with 



E[Cl] 



-h-n 



(41) 



(42) 



the energy of the spin in the external field. It turns 
out that Eq. (|4"Tj) holds for any effective field h = 
—dE[ti]/dCl, and in particular for the case that E[Cl] is 
quadratic in the components of Cl as is often used to 
model magnetic anisotropy. 

It is important to realize that the equilibrium prob- 
ability distribution has precisely this form because of 
the fluctuation-dissipation theorem, which ensures that 
dissipation and fluctuations cooperate to achieve ther- 
mal equilibriu m 13 ^ 4 . Finally, it should be noted that 
this derivation of the stochastic Landau-Lifschitz-Gilbert 
equation from a microscopic starting point circumvents 
concerns regarding the phenomenological form of damp- 
ing and fluctuation-dissipation theorem, which is subject 
of considerable debat o 22 ' 23 . 



IV. NONZERO BIAS VOLTAGE 

In this section we consider the situation that the chem- 
ical potential of the left lead is given by /Ul = \i + \e\V, 
with \e\V > the bias voltage in units of energy, and 
fi = (Ur, the chemical potential of the right lead. Using 
the general expressions given for the response functions 
derived in Sec. HTJ it is easy to see that the imaginary 
part of the retarded and advanced components of the 
response functions are no longer related to the Keldysh 
component by means of the fluctuation-dissipation theo- 
rem in Eq. (|3"4"]) . See also the work by Mitra and Millis 28 
for a discussion of this point. As in the previous section, 
we proceed to determine the low-frequency behavior of 
the response functions. 

Using Eqs. ([26]), (|27]), and (J30J) we find that the re- 
tarded and advanced components of the response func- 
tion are given by 

n(±) (<0 = [A 2 0" + W) + A 2 (/j)] e . (43) 



In this expression we have omitted the energy- 
independent part and the contribution following from the 
principal-value part of the energy integral because, as we 
have seen previously, these do not contribute to the final 
equation of motion for the direction of the spin. Follow- 
ing the same steps as in the previous section, we find 
that the damping kernel in the general Langevin equa- 
tion in Eq. (|25p reduces to a Gilbert damping term with 
a voltage-dependent damping parameter given by 



*M = -^[A 2 ^+\e\V)+A 3 0x)] 



cv 



1 + 



\e\V 
A* 



(44) 



This result is physically understood by noting that the 
Gilbert damping is determined by the dissipative part of 
the response function I[( + )(e). In this simple model, this 
dissipative part gets contributions from processes that 
correspond to an electron leaving or entering the system, 
to or from the leads, respectively. The dissipative part 
is in general proportional to the density of states at the 
Fermi energy. Since the Fermi energy of left and right 
lead is equal to fj, + \e\V and /i, respectively, the Gilbert 
damping has two respective contributions corresponding 
to the two terms in Eq. (j44|) . 

Note that the result that the Gilbert damping param- 
eter initially varies linearly with the voltage is in agree- 
ment with the results of Katsura ef al£l, although these 
authors consider a slightly different model. In the inset of 
Fig. [T] we show the Gilbert damping parameter as a func- 
tion of voltage. The parameters taken are A/ J = 0.1, 
J c = J, (i/J = 1 and n/(k B T) = 100. 

Although we can no longer make use of the fluctuation- 
dissipation theorem, we are nevertheless able to deter- 
mine the fluctuations by calculating the low-energy be- 
havior of the Keldysh component of the response function 
in the nonzero-voltage situation. It is given by 



n K (e) 



dc' 



-A 2 (e'){[N( PL -e') + N( m -e')} 

(2?r) 

x [N(e'-fx L )+N(e'- m )]} . 



We define an effective temperature by means of 

_ zII K (e)A 2 



k B T cS (T, V) 



2a(V) 



(45) 



(46) 



This definition is motivated by the fact that, as we 
mention below, the spin direction obeys the stochas- 
tic Landau-Lifschitz-Gilbert equation with voltage- 
dependent damping and fluctuations characterized by the 
above effective temperature* 3 -. From the expression for 
a(V) and II K (e) we see that in the limit of zero bias 
voltage we recover the equilibrium result T c « = T . In 
the situation that \e\V is substantially larger than k B T, 
which is usually approached in experiments, we have that 



k B T cfi (T, V) 



\e\V , k B T 



(47) 
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which in the limit that \e\V ^> k-gT becomes indepen- 
dent of the actual temperature of the leads. In Fig. Q] 
the effective temperature as a function of bias voltage is 
shown, using the expression for II K (e) given in Eq. (|45[) . 
The parameters are the same as before, i.e., A/J = 0.1, 
Jc = J, n/J = 1 and li/(knT) = 100. Clearly the ef- 
fective temperature changes from T e e = T at zero bias 
voltage to the asymptotic expression in Eq. (|47[) shown by 
the dashed line in Fig. [T] The crossover between actual 
temperature and voltage as a measure for the fluctuations 
is reminiscent of the theory of shot noise in mesoscopic 
conductors 44 . This is not surprising, since in the single- 
site model we use the noise in the equation of motion ul- 
timately arises because of fluctuations in the number of 
electrons in the single-site system, and is therefore closely 
related to shot noise in the current through the system. 
Foros et a/.— calculate the magnetization noise arising 
from spin current shot noise in the limit that \e\V 3> k^T 
and |e|V <C k^T. In these limits our results are similar 
to theirs. 

With the above definition of the effective temperature 
we find that in the nonzero bias voltage situation the spin 
direction obeys the stochastic Landau-Lifschitz-Gilbert 
equation, identical in form to the equilibrium case in 
Eqs. and P0"|) , with the Gilbert damping parame- 
ter and temperature replaced according to 

a a -> a(V) ; 

T -> T eG (T,V) . (48) 

Moreover, the transport-steady-state probability distri- 
bution for the direction of the spinis a Boltzmann distri- 
bution with the effective temperature characterizing the 
fluctuations. 



V. DISCUSSION AND CONCLUSIONS 

We have presented a microscopic derivation of the 
stochastic Landau-Lifschitz-Gilbert equation for a semi- 
classical single spin under bias. We found that the 
Gilbert damping parameter is voltage dependent and to 
lowest order acquires a correction linear in the bias volt- 
age, in agreement with a previous study for a slightly 
different model^l. In addition, we have calculated the 
strength of the fluctuations directly without using the 
fluctuation-dissipation theorem and found that, in the 
low-frequency regime, the fluctuations are characterized 
by a voltage and temperature dependent effective tem- 
perature. 

To arrive at these results we have performed a low 
frequency expansion of the various correlation functions 
that enter the theory. Such an approximation is valid as 
long as the dynamics is much slower than the times set 
by the other energy scales in the system such as temper- 
ature and the Fermi energy. Moreover, in order for the 
leads to remain in equilibrium as the spin changes direc- 
tion, the processes in the leads that lead to equilibration 



have to be much faster than the precession period of the 
magnetizationspin. Both these criteria are satisfied in 
experiments with magnetic materials. In principle how- 
ever, the full Langevin equation derived in Sec.[TT]also de- 
scribes dynamics beyond this low-frequency approxima- 
tion. The introduction of the effective temperature relies 
on the low-frequency approximation though, and for ar- 
bitrary frequencies such a temperature can no longer be 
uniquely defined^ 8 -. 

An effective temperature for magnetization dynam- 
ics has been introduced before on phenomenological 
ground o 36 i 37 ' 38 . Interestingly, the phenomenological ex- 
pression of Urazhdin et al£^, found by experimentally 
studying thermal activation of current-driven magneti- 
zation reversal in magnetic trilayers, has the same form 
as our expression for the effective temperature in the 
large bias-voltage limit [Eq. (|4"T]) ] that we derived micro- 
scopically. Zhang and Li 3 -^, and Apalkov and Visscher 38 , 
have, on phenomenological grounds, also introduced an 
effective temperature to study thermally-assisted spin- 
transfer-torque-induced magnetization switching. In 
their formulation, however, the effective temperature is 
proportional to the real temperature because the current 
effectively modifies the energy barrier for magnetization 
reversal. 

Foros et al£^ consider spin current shot noise in the 
large bias- voltage limit and find for sufficiently large volt- 
age that the magnetization noise is dominated by shot 
noise. Moreover, they also consider the low bias-voltage 
limit and predict a crossover for thermal to shot-noise 
dominated magnetization fluctuations. Our main result 
in Fig. Q] provides an explicit example of this crossover 
for a simple model system obtained by methods that are 
easily generalized to more complicated models. In the 
experiments of Krivorotov et the temperature de- 
pendence of the dwell time of parallel and anti-parallel 
states of a current-driven spin valve was measured. At 
low temperatures k^T < \e\V the dwell times are no 
longer well-described by a constant temperature, which 
could be a signature of the crossover from thermal noise 
to spin current shot noise. However, Krivorotov et al. 
interpret this effect as due to ohmic heating, which is 
not taken into account in the model presented in this 
paper, nor in the work by Foros et al£^. Moreover, in 
realistic materials phonons provide an additional heat 
bath for the magnetization, with an effective tempera- 
ture that may depend in a completely different manner 
on the bias voltage than the electron heat-bath effec- 
tive temperature. Nonetheless, we believe that spin cur- 
rent shot noise may be observable in future experiments 
and that it may become important for applications as 
technological progress enables further miniaturization of 
magnetic materials. Moreover, the formalism presented 
here is an important step in understanding magnetization 
noise from a microscopic viewpoint as its generalization 
to more complicated models is in principle straightfor- 
ward. Possible interesting generalizations include mak- 
ing one of the leads ferromagnetic (see also Ref. [4q). 
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Since spin transfer torques will occur on the single spin 
as a spin-polarized current from the lead interacts with 
the single-spin system, the resulting model would be a 
toy model for microscopically studying the attenuation 
of spin transfer torques and current-driven magnetiza- 
tion reversal by shot noise. Another simple and use- 
ful generalization would be enlarging the system to in- 
clude more than one spin. The formalism presented here 
would allow for a straightforward microscopic calcula- 
tion of Gilbert damping and adiabatic and nonadiabatic 
spin transfer torques which are currently attracting a lot 
of interest in the context of current-driven domain wall 
motionSiLSi i 9 ! 10 ! 11 ! 12 . The application of our theory in the 



present paper is, in addition to its intrinsic physical inter- 
est, chosen mainly because of the feasibility of analytical 
results. The applications mentioned above are more com- 
plicated and analytical results may be no longer obtain- 
able. In conclusion, we reserve extensions of the theory 
presented here for future work. 
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